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Abstract 

We review recent results about the derivation and the analysis of two Hartree-Fock- 
type models for the polarization of vacuum. We pay particular attention to the variational 
construction of a self-consistent polarized vacuum, and to the physical agreement between 
our non-perturbative construction and the perturbative description provided by Quantum 
P-t Electrodynamics. 

^ ■ 1 Introduction 

During the past century, the classical picture of the vacuum as an empty object was challenged by 
a series of theoretical advances and experimental observations including the measurement of the 
Lamb shift |29] and the derivation of the Casimir effect |3] (see also jl]). The description of the 
vacuum as a complicated fluctuating system emerged. An intuitive picture for this system can be 
derived from the observation by Blackett and Occhialini [2] of the creation of electron-positron 
pairs when one provides a sufficient amount of energy to the vacuum. Since the energy-time 
uncertainty principle |36) allows important fluctuations of energy during short time intervals, 
nothing prevents the vacuum from being the place of permanent creations and annihilations of 
pairs of virtual particles. 

This phenomenon affects the interactions between physical particles. It is in particular at the 
origin of the vacuum polarization. In presence of an external electromagnetic field, the virtual 
electron-positron pairs modify the distributions of charges and currents originally generated by 
| the fields. The corrections with respect to the original distributions are computed using Quantum 

Electrodynamics (see e.g. [40, 19, 10J). This perturbative theory provides their value in terms 
of a power series with respect to the Sommerfeld fine-structure constant 
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X: 

Unfortunately, divergences appear at any order in the computations. In order to obtain a well- 
defined value, one has to appeal to complicated renormalization procedures, which provide cor- 
rections in an extremely accurate agreement with physical experiments in spite of their intricacy 
(see e.g. [15J). 
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Our goal in this survey is to describe recent mathematical results concerning simplified models 
for the vacuum polarization. The main difference with Quantum Electrodynamics lies in the 
non-perturbative nature of the models. The fine-structure constant a being fixed, the polarized 
vacuum is constructed by variational arguments. The main difficulty lies in the choice of the 
approximations to make in order to allow such a construction and to guarantee its relevance with 
respect to the perturbative computations of Quantum Electrodynamics. 

The results in this survey are based on a seminal paper by Chaix and Iracane [5] (see also 
[6]), which provides a mean-field framework for the analysis of the vacuum polarization. In this 
setting, the system under consideration is composed of the physical electrons and positrons, cou- 
pled with virtual ones which give account of the polarized vacuum. All of them interact with 
photons, the interactions being instantaneous. The electrons and positrons are described using 
an Hartree-Fock approximation. Under this approximation, and using renormalization proce- 
dures which are quite standard in the context of Quantum Electrodynamics, it becomes possible 
to define rigorously an energy for the system. The polarized vacuum is then mathematically 
constructed as a critical point of this energy. One can compare the expansions of its charge and 
current densities with respect to the fine-structure constant a to the ones provided by Quantum 
Electrodynamics, and check, at least in the purely electrostatic case, that the approximations 
made in order to fashion the simplified models are quite reasonable. 

This is in brief what we are going to describe in this survey. In the first section, we focus 
on the derivation of the two Hartree-Fock models that we are going to analyze mathematically. 
We provide some further details on the notion of polarized vacuum, as well as on the nature 
of an Hartree-Fock approximation. The second section is devoted to the mathematical analysis 
of our first model, the so-called Bogoliubov-Dirac-Fock model, which corresponds to the purely 
electrostatic case. We give rigorous existence results for the polarized vacuum, and check (in 
the reduced case) the relevance of our definitions with respect to Quantum Electrodynamics 
by computing an asymptotic expansion of the total charge density with respect to the fine- 
structure constant a. In the third section, we consider the Pauli-Villars regulated model which 
takes account of the electromagnetic interactions with the photons. Our main results about 
this further model deal with its rigorous definition and with the mathematical construction of a 
polarized vacuum. 

2 Derivation of the Hartree-Fock models for the vacuum polar- 
ization 

2.1 The picture of the Dirac sea 

In relativistic quantum mechanics, individual electrons are represented by spinors tp £ L 2 (R 3 , 
C 4 ). This description originates in the works of Dirac [7J |8, 9J, who introduced the formula 

E c (lJ)) = Pm,o(^),'0)L 2 (M3 ) C 4 ), 

for the computation of the kinetic energy of a relativistic electron. The free Dirac operator D m fi 
is defined as 

D m ,o = to • (— iV) + mc 2 (3, 

where K, c and m stand respectively for the reduced Planck constant, the speed of light (in the 
free vacuum), and the (bare) mass of an electron. Without loss of generality, we can make a 
choice of units such that h = 1 and c = 1. In the sequel, we adopt this choice, so that we drop 
the dependence on h and c of the operator D m $. 
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The Dirac operator is a self-adjoint operator on L 2 (M 3 , C 4 ), with domain ff 1 (M 3 , C 4 ) (see e.g. 
|48j ) . The Dirac matrices a = (0:1,0:2,0:3) and (3 are given by the formulae 

ak= {a k ^ and /3= (o -I 2 )' 
where the Pauli matrices o"i, a 2 and <x 3 are equal to 

The Dirac matrices are designed so as to satisfy the identity 

D 2 mfi = -A + m 2 I. 

As a result, the spectrum of the free Dirac operator splits into two intervals according to the 
expression 

c(An,o) = (-OO, -TO] U [TO, +00). 

The spectrum is interpreted as the possible levels of kinetic energy for a relativistic electron, so 
that nothing prevents the kinetic energies to be arbitrarily negative. Such a phenomenon has 
never been observed in practice. 

Dirac [7J |8j |9] by-passed the difficulty introducing the picture of the Dirac sea. In the free 
vacuum, an infinite number of virtual electrons completely fill in the levels of negative energy. The 
possible levels of energy for "physical" electrons are positive. The picture amounts to claiming 
that the free vacuum is not represented by a vanishing quantity, but instead is identified to the 
negative spectrum of the free Dirac operator. In the following where we will consider electrons 
in an Hartree-Fock state (see Subsection 12. 2\ , the free vacuum will be more precisely identified 
to the negative spectral projector 

P m,0 = X(-oo,0](An,o)- 
In presence of an additional external electromagnetic four-potential 

^4-ext = (^extj ^4-ext)) 

the nature of the vacuum is modified by the interactions between the virtual electrons and 
the external field. In the simplified Furry picture |14j . the vacuum is described through the 
introduction of the electromagnetic Dirac operator 

D m ,eA^ t = ol ■ (-z'V - eA ext ) + to/3 + eKxt, 

which is written here in physical units such that the vacuum permittivity £q is equal to 1/47T. 
The so-called dressed vacuum is identified to the negative spectral projector 

P m,A c ^t = X(-oo,0](An,eA oxt )- 

The virtual electrons of the dressed vacuum span the range of the projector P^™a ■ They 
have no reason to remain identical to the virtual electrons in the free (or bare) vacuum which 
span the range of the projector -P~ - ^ n g enerai ) the charge density of the vacuum does not 
remain constant. The dressed vacuum is polarized. In the case where the external field is strong 
enough, a positive eigenvalue can even appear in the spectrum of the electromagnetic Dirac 
operator D mt A cxt , creating an hole in its negative spectrum. A (physical) electron-positron pair 
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is produced, the electron being identified to the positive eigenvalue, while the hole is identified 
to the positron. 

The Furry picture provides a good approximation for the polarization of the vacuum when the 
external fields are not too strong. In practice, the non-constant charge density of the polarized 
vacuum modifies the electromagnetic field. The virtual electrons react to the corrected field 
which in turn affects the nature of the polarization. In order to describe it more accurately, one 
has to look for more sophisticated models. 



2.2 The Hartree-Fock approximation 

In the Physics literature, vacuum polarization is described using Quantum Electrodynamics 
which provides extremely accurate computations for the charge and current densities of the 
polarized vacuum. On the mathematical level, this theory is far from being completely under- 
stood, in particular, due to the perturbative nature of its computations. In this subsection, we 
present a set of approximations which make possible the construction of non-perturbative models 
for the vacuum polarization. The main one consists in describing the electronic structures as 
Hartree-Fock states. 

The simplest way to introduce the Hartree-Fock approximation is probably to come back to 
the description of N classical electrons around a positive density of charge v. In this situation, 
the electronic structure is described through the iV-body Hamiltonian 



N 

tN 
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which acts on the space L^((M 3 ) Ar , C) of electronic states Vl/ which are anti-symmetric with respect 
to the permutations of the variables Xi (in order to guarantee the validity of the Pauli exclusion 
principle), and with a density \^>\ 2 which depends symmetrically on the variables Xj (so as to 
handle with undistinguishable electrons). The possible electronic structures correspond to the 
eigenfunctions of the Hamiltonian , the ground state structure corresponding to the minimal 
eigenvalue. The analysis of H„ is rather involved (see |28| I3"T] for more details), so that several 
approximations have been suggested to simplify the description. 

The Hartree-Fock approximation |25[ [T3] restricts the analysis to the Hartree-Fock states '3/ 
which write as Slater determinants 

■$(xi,...,x N ) = ttiA---A¥;v(xi,...,a:jv) = -= det (*i(arj))j 



/]y[ \ = l ^'n<i,j<N> 



where ^i, are N orthogonal wavefunctions in L (M , C). The Hartree-Fock states 

A • • • A \&at are the less correlated electronic structures. Their energy is given by the expression 

N / r \ „2 /• /• u, ».M2 



2m ~{ V Jr 3 ) 2 jr? Jr 3 \ x - y\ 

7 / 

Jr s Jr 



u{x)My) dxdy+ e 4 [ [ Mx)My) dxdy. 



i=l 

//j " ' 2 ./;::; ./;::; \.V - 1/ 

In this formula, 7^ refers to the one-body density operator with integral kernel 

N 



jqt(x,y) =J2^i(x)^i(y). 



i=l 
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In other words, 7^ is the orthogonal projector on the linear space spanned by the wavefunctions 
^j. Concerning the charge density p^, it is equal to 



N 



i=l 

In particular, the Hartree-Fock energy only depends on the one-body density operator 7^ through 
the identity 

E u (7*) =— tr(-A7^) - — / / — —dxdy 

, , , . ( 2 ) 



2?n 2 J R 3 J R 3 \x-y 

2 



j , f 'M£lM dxdy + & - f f P ^ x)p ^ y) dxdy, 
Jus J R 3 \x-y\ 2 J R3 J R 3 \x - y\ 

where tr(— A 7^) is the trace of the finite rank operator —A 7^. As a consequence, the Hartree- 
Fock electronic structure can be computed in terms of an orthogonal projector, the Hartree-Fock 
ground state being identified to the projector which minimizes the energy E^ F among all the 
possible projectors. 

The picture is similar in the relativistic case. In Quantum Electrodynamics, the (formal) 
Hamiltonian describing an electronic structure in the presence of a (classical) external electro- 
magnetic four-potential A cxt = (Koctj^ext) may be written in the Coulomb gauge as 



= J **(x) a ■ (-iV - eA(x) - eA cxt (x)) + m/3 

+e f V cxt (x)p(x)dx + — I I pM^M dxdy + — [ I curl A(x)\ 2 dx. 
J 2 J J \x-y\ 8tt J 



Hf(x) dx 

2 , , , , (3) 



In this expression, ^f(x) refers to the second quantized field operator which annihilates an electron 
at x, while the vector A(x) is the magnetic field operator for the photons. The density operator 
p{x) is defined as 

p(x)= 1 -Y / [K(x)^Ax)], (4) 

CT=1 

where [a, b] = ab - ba (see [26l [27J S51 S31 [23l [32] for more details about the Hamiltonian M Acxt ). 

The Hamiltonian EI" 4 ™' acts on the Fock space T = T c where F e is the fermionic Fock 

space for the electrons and .Fph is the bosonic Fock space for the photons. As in the usual case, 
the main approximation in order to derive Hartree-Fock models for the description of relativistic 
electrons consists in restricting the Hamiltonian to states of the special form 

f2 = ^HF ® ^Coh, 

where Ohf is an electronic Hartree-Fock state characterized by its one-body density matrix 

j(x,y) = (**(x)*(l/))nap. 
and l^coh is a coherent state characterized by its classical magnetic potential 

A(x) = {A(x))n Cob - 

In practice, this amounts to representing the electrons by an operator 7 as in classical Hartree- 
Fock theory, while the photons are described by a magnetic potential A, which is nothing more 
than a classical vector field. The main difference with the classical case lies in the fact that 
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the operator 7 does not only take into account the physical electrons, but also the virtual ones 
representing the vacuum. In general, it is not anymore an orthogonal projector with finite rank, 
but with infinite one, which causes the apparition of divergences in the computation of the energy. 

Up to a universal constant which diverges in infinite volume, this energy is equal to 

«^'(7,A)= t r( Dra , e(w+ ,„, ) ( 7 -l/2)) + f / / 

v / z j K 3 j^3 \x — y\ 

/ / ^ — . 7 dgdy + — / | curl ,4(x)| 2 dec. 

2 Jr3 y R a |sc - y| 8vr J R3 

This expression is similar to the one for the classical Hartree-Fock energy in ([2j. The kinetic 
energy is now equal to tr(D mj o(7 — 1/2)), which amounts to replacing the classical operator 
—A by its relativistic version D m fi. One also recovers the presence of the so-called direct and 
exchange electrostatic terms in the first, respectively second, line. The main difference lies in 
the property that the energy is not directly written in terms of the density operator 7, but in 
terms of the difference 7 — 1/2. This is a consequence of the charge-conjugation invariant choice 
(j3J for the density operator p{x) (see [23] for more details). This property is particularly helpful 
when one attempts to give a rigorous meaning to the expression in ([5]). 

The energy in (|5]) presents major divergences. Since the operator D m $ is unbounded and 
7 — 1/2 has infinite rank when 7 is a projector, the kinetic energy is not well-defined. Moreover, 
since 7 — 1/2 is not a compact operator, defining its kernel (7 — l/2)(x,y), and thereafter the 
density p 7 _i/2> is a further challenge. 

Similar divergences appear in Quantum Electrodynamics, in which regularization techniques 
have been developed in order to by-pass the difficulty These techniques mainly rely on the 
introduction of an ultraviolet cut-off A in the model. In the sequel, we describe how it is 
possible to adapt two regularization techniques to provide a rigorous meaning to the energy in 
(0) and construct a consistent Hartree-Fock model for the vacuum polarization. This leads to the 
construction of two different models. The first one is the Bogoliubov-Dirac-Fock model which 
does not take into account the effects due to the presence of photons, as well as of a possible 
external magnetic field. In the second one, handling with these effects is in contrast possible due 
to the use of a more ingenious regularization, the so-called Pauli-Villars regularization. 



3 The Bogoliubov-Dirac-Fock model 
3.1 Derivation and functional framework 

The Bogoliubov-Dirac-Fock model was introduced by Chaix and Iracane in [5] (see also [6]). 
It only takes into account kinetic and electrostatic aspects. The external magnetic potential 
A ext and the magnetic potential A for the photons are set equal to 0. Concerning the external 
electrostatic potential V ex t, it is induced by an external charge density v according to the Coulomb 
formula 

V cxt (x) = -e [ r^-dy. 



/R3 \x - y\ 

As a result, the formal Hartree-Fock energy for the one-body density operator 7 reduces to the 
expression 



£gp(7) =tr (D m)0 ( 7 - 1/2)) -all h l/^M dxdy 



2 7 7 \x — y\ 2 J J \x — y\ 
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where we have introduced the (bare) Sommerfeld fine-structure constant a = e 2 . Recall that in 
this formula, the function (7 — l/2)(x,y) refers to the (formal) kernel of the operator 7 — 1/2, 
while the density p 7 _i/2 is (also formally) defined as 

Pj-i/2(x) = (7 - l/2)(x,x). 

One can introduce a reduced version of the energy omitting the exchange electrostatic term 
according to the formula 

fr V(7) = tr - 1/2)) -aff ^M^M dldy + °[f ^ 

J J \x -y\ 2 J J \x-y\ 

(7) 

In the sequel, we restrict our attention to the reduced energy in order to simplify the analysis. 
We will point out the results which remain available for the energy with exchange term (see 
|23[ [T2| [32] and references therein for more details). 

The critical points for the energy are solutions to the self-consistent equations 

[ 7 ,D 7 ] =0, 

where the Fock operator is defined as 

£> 7 = D m>0 + a(p y _ 1/2 - v) * -r-r ■ 

An orthogonal projector 7* which minimizes the energy among all the possible projectors, 
is solution to the equation 

7* = Xi-ocfi^D-y,)- (8) 

The projector 7* is interpreted as the polarized vacuum according to the picture of the Dirac sea. 
The spinors of the virtual electrons in the polarized vacuum span the range of 7*. In view of the 
self-consistent equation (|8]), they completely fill in the negative spectrum of the Fock operator 

The model is fashioned to allow for a description of the electronic structure with iV physical 
electrons around the charge density v. When 7 is an orthogonal projector with finite rank, its 
rank, or alternatively its trace, provides the number, or alternatively the total charge, of the 
electrons in the structure represented by 7. The ground state structure is defined as a minimizer 
7at of the energy £^ F in the charge sector with charge N, that is among the orthogonal projectors 
with trace equal to N. Such a minimizer (formally) satisfies the self-consistent equation 

IN = X(-oo,n N ]{D 1N ), (9) 

where /ijy is the Fermi level of the electronic structure described by the operator 7jy. Following 
again the picture of the Dirac sea, the projector 7™° = X(-oo.o](-^7jv) describes the virtual 
electrons of the vacuum polarized by the charge density v and the N physical electrons. They 
are themselves given by the orthogonal projector 7^ = X(o,[i N ]{D^ N )- In particular, their energy 
is positive. 

At this stage, it is necessary to emphasize that most of the definitions in the previous discussion 
are only formal. For example, any orthogonal projectors of finite trace must have a finite rank 
which cannot be consistent with the identity in (j9j). However, it is possible to provide a rigorous 
meaning to the reduced Bogoliubov-Dirac-Fock model (defining for instance a notion of trace for 
projectors with infinite rank) by invoking regularization techniques. 

In this direction, Hainzl, Lewin and Solovej |24] noticed that the model is well-defined in a 
finite-dimensional setting. They suggested to restrict the analysis to a box Cl = [— L/2,L/2[ 3 
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of size L, with periodic boundary conditions, and to introduce an ultraviolet cut-off A in the 
Fourier domain. This amounts to assuming that the operators 7 act on the finite-dimensional 

space 

where 1Zl,A = {k £ (2-k/L) Z 3 , s.t. \k\ < A}. In this periodic framework, the reduced Hartree- 
Fock energy takes the form 

ZpcAl) = tr (^,0(7 - V 2 )) - a / / P-y~i/2(x)u L (y)G L (x - y) dxdy 

Jc L Jc L , , 

a f f { ' 

+ 2 / Pi-i/2(, x )P-y-i/2(y)G L (x -y)dxdy. 

In this definition, the function vl is a periodic version of the charge density v. It is equal to 

^ 3 



///(./•) -_[-L^!L\ ^ V(k)expi(k,x) 



where v refers to the Fourier transform of v. In order to get a chance to define properly the 
Coulomb potentials in (|10p . it is natural to assume that u belongs to the Coulomb space 



C(M 3 ) = {/ € Ll c (M 3 ,C), s.t. \\f\lt = M- dk < +00}. 
The periodic Coulomb kernel is similarly given by 

G L (z) = ^ ( |Zj2 exp«'(A;,^) M 3 J +K L , 

fcG(27r/L) Z3\{0} ^' 1 ' 

where the constant Kl is fixed so that Gl is positive. Finally, the charge density Pj-1/2 is 
defined as 

Pi-i/2( x ) = tr c 4 (7 - l/2)(x,x) = -^3 ^ tr C 4 [7 - l/2](j,k) expi(j - k,x) R3 , 

where (7 — l/2)(x,y) is the kernel of the (finite-rank) operator 7 — 1/2. 

In this setting, the periodic reduced Hartree-Fock energy is well-defined on the convex hull of 
orthogonal projectors on Hl,A, which is defined as 

G L ,A = {76 C(Hl,a), s.t. 7* = 7 and < 7 < /}. 

Moreover, the energy owns a minimizer on Gl,A, which can be interpreted as the polarized 
vacuum 

Proposition (|24]). Let m, a, A and L be positive numbers and consider a function v £ C(R 3 ) 
such that v is continuous on R 3 . There exists a minimizer 7^ to the minimization problem 

mf e;M- 

When v = 0, the minimizer 7^ is the negative spectral projector P~ L °f ^ e restriction of the 
free Dirac operator D m ^ to the finite- dimensional space Hl,a- 



1 Extending the minimization problem under consideration to a convex hull like Gl,a is standard in Hartree- 
Fock theory. The construction of minimizers is simplified, and it turns out that they still solve the initial problem 
(see e.g. [HI]). 
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The operator 7^ is identified to the polarized vacuum corresponding to the previous periodic 
setting. A simple way to describe the polarized vacuum in the full space is to consider the 
thermodynamic limit L — > +00 of the minimizers 7^. It turns out that the limit exists (when the 
ultraviolet cut-off A is fixed) and that it can be described using the reduced Bogoliubov-Dirac- 
Fock energy (see [24J). 

This energy is defined according to the formula 

%dp(Q)=V (D m , Q)-af v -fE9f dxdy + ^f f E9ME9M dxdy . (11) 

m '° 7r3 \x-y\ ZJrzJrz \x-y\ 

The difference with respect to the reduced Hartree-Fock energy in ([7J lies in the choice of a 
reference projector to compute the energy. More precisely the energy is not directly expressed 
in terms of the projector 7, but in terms of its difference 

with the negative spectral projector P~ of the free Dirac operator. In other words, the free 
vacuum energy is set equal to and all the energies are computed with respect to this reference. 

Hainzl, Lewin and Sere |20] constructed a functional framework in which the reduced Bogoliu- 
bov-Dirac-Fock energy is well-defined and bounded from below. They observed the necessity to 
conserve an ultraviolet cut-off A in the model (see |21] ) assuming that the operators Q do not 
act on L 2 (R 3 ,C 4 ), but on the space 

U K = {* € L 2 (R 3 ,C 4 ), s.t. supp(§) C 5(0, A) j. 

This assumption presents the major advantage to transform the operator D m $ in a bounded 
operator. We will see in Subsection 13.31 below how to eliminate the ultraviolet cut-off A from 
this description of the polarized vacuum. 

More precisely, the operators Q belong to the function space 

6 A = {TG 6 2 (H A ), s.t. (P- TP- , (/ - P- )T(/ - P- )) € &i(n A ) 2 }, (12) 

where &i(T~L\) and &2('Ha) are the spaces of trace-class, respectively Hilbert-Schmidt, operators 
on %a- The trace-class conditions in (|12p originate in the property that a formal minimizer 
of the energy S^BDF * s m g enera l n °t trace class (see Theorem [2] below) . As a consequence, it 
is not so easy to define properly its kinetic energy. This is done extending the usual definition 
through the formula 

tr ^,o ( T ) = tr ( P m,o TP m,o) + * ((/ - P-, )T(I ~ P„, )) , (13) 
which is well-defined when T is in ©9 . Notice here that the quantity tr p - (T) is equal to the 

m,0 

trace of T when T is trace class. 

Concerning the charge density pq, its definition relies on the introduction of the ultraviolet 
cut-off A. Since the operator Q is Hilbert-Schmidt, it owns a kernel Q(x,y), which belongs 
to the space L 2 (M 3 x M 3 ). Moreover, its Fourier transform is supported into the product set 
B(0, A) x 5(0, A) due to the presence of the cut-off. In particular, it is a smooth function so 
that the charge density pq can be defined according to the identity 

p Q (x) = tr C 4 Q(x,x). 

With this last definition at hand, one can define properly the reduced Bogoliubov-Dirac-Fock 
energy H. 

2 The original Bogoliubov-Dirac-Fock energy with exchange term is also well-defined on a set similar to (3a- 
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Proposition ([20"]). Let A > 0, a > and v £ C(R 3 ). The energy S^bdf ^ s well-defined on the 
set 6°. 

The introduction of the reduced Bogoliubov-Dirac-Fock model is justified by taking the ther- 
modynamic limit L — > +oo in the periodic model introduced by Hainzl, Lewin and Solovej in 
[24] . As a matter of fact, they established the following theorem (which remains available with 
a few minor changes in the case of the Bogoliubov-Dirac-Fock model with exchange term). 

Theorem ([24]). Let m, a and A be positive numbers and consider a function v G C(R 3 ) such 
that v is continuous on R 3 . 
(i) When v = 0, we have 

In particular, 
and 



ll = P m ,0,L- 

H-P^o\\ c{nLA) ^0, (14) 

£per(7£) -> +00, (15) 

as L —> +oo. 

(ii) There exists an operator in & A such that, up to a possible subsequence, we have 

Ql( x ,v) = iU x ,y) -il( x iy) -> Q*( x ,y), (16) 

as I — > +oc, uniformly on any compact subset o/R 3 x R 3 . 

(iii) The operator is a minimizer of the reduced Bogoliubov-Dirac-Fock energy £^ BDF on the 
set 

Q A = [Q G e% s.t. Q* = Q and - P m0 <Q<I- P mfi ). 
Moreover, we have 

%M " £per(7i) -> ^BDf(Q*) = min {KbDp(Q), ^ Q A }. (17) 

as L — )• +oo. 

The choice of the projector P~ as reference for a model in the full space is justified by the 
convergences in (|14p and (|15p . In the thermodynamic limit L — > +oo, the free vacuum is given 
by P^o according to the original picture of the Dirac sea. Moreover, one has to subtract its 
infinite energy in order to handle with a reasonable model. 

In this case, the convergences in (|16p and (|17p show that the polarized vacuum is actually 
described by the reduced Bogoliubov-Dirac-Fock model. More precisely, the polarized vacuum is 
identified to the projector 7* = P~ + Q* corresponding to the minimizer Q* of the energy £^ BDF 
on the set Q\. Let us recall now some elements about the construction of such a minimizer, as 
well as its main properties. 

3.2 Construction of the polarized vacuum and the electronic structures with 
electrons 

The existence of the minimizer Q* was established by Hainzl, Lewin and Sere in |20[ |21] . 
Theorem ([HJEI]). Let A > ; a > and v G C(R 3 ). 

{£) The energy £^ BDF is bounded from below on the set Q\. More precisely, we have 

3bdf«) + ?/ / V -P^-dxdy>Q, (18) 
2 Jr3 J R 3 \x - y\ 
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for any Q £ Q\. When u = 0, inequality (jl8|) terns into an equality if and only if Q = 0. 

energy £^QY)p owns a minimizer on Q\, which is solution to the self- consistent 

equation 

7* = Q* + Pm,o = X(-oo,o)( D *) + 8*, (19) 
where the Fock operator D* is equal to 

D* = D mfi + a(pQ„ - v) * r^r, 

u;/ii/e J* refers to a finite-rank operator with range in the kernel of . 

[Hi) The charge density pQ t is uniquely determined. When the external charge density v satisfies 
the condition 

11 1 1 
2 a 7T6Q!||z/||c < m ! , (20) 

the minimizer is also unique. In this case, the operator 5* is equal to 0, and 

N* = tip- (Q*) = 0. (21) 

Remark. The previous theorem extends with a few minor changes to the Bogoliubov-Dirac-Fock 
model with exchange term (see |20[ 124] ). In this case, inequality (|18p was first proved in [lj. 

Once again, the fact that Q = is the unique minimizer of the energy <£^bdf ^ s consistent 
with the picture of the Dirac sea for the free vacuum. In presence of an external charge density 
v, the polarized vacuum is described by the operator 7* = + P~ , which is not necessarily an 
orthogonal projector due to the presence of the finite-rank operator 5*. This defect is a common 
drawback of the reduced Hartree-Fock models (see e.g. [46j). It is possible to withdraw the 
operator 5* from the self-consistent equation satisfied by a minimizer of the Bogoliubov-Dirac- 
Fock energy with exchange term (see |20| I24|). 

Under condition (|20p . the kernel of the Fock operator is necessarily equal to 0, so that the 
operator 5* identically vanishes. In this case, one recovers the picture of the Dirac sea in which 
the polarized vacuum is identified to the negative spectral projector 7* of the Fock operator D*. 
It follows from (|2ip that the polarized vacuum is globally neutral. The trace of Q* is indeed 
interpreted as the charge of the electronic structure represented by 7* since, at least formally, 
the charge of the negative projector P~ , which gives account of the free vacuum, must be equal 
to 0. This is consistent with the observation that a weak electrostatic potential cannot produce 
physical electrons in the vacuum. 

In another direction, one can ask for the behaviour of the minimizer Q* when A — > +00. It 
turns out that the model collapses in this limit. 

Proposition (|21]). Let a > and v G C(M 3 ). We denote by a minimizer of the energy 
^rBDF 071 2a f or a P os iti ve number A. Then, 

\\\D m ,o\^Qt ||g 2 and a\\p Q A - v\\ c 0, 
as A — > +00. In particular, 

/ (22) 
2 Jr3 Jrs \x - y\ 

as A — > +00. 
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In the limit A — > +00, the operators converge to 0, while their charge densities Pqa tend 
to the function v. The limit density is not the charge density of the limit operator. In other 
words, the model does not remain physically consistent in the limit A — > +00. In particular, 
taking this limit is not a way to eliminate the ultraviolet cut-off A. This property is related to 
the Landau pole phenomenon which was originally described by Landau [30j, and Landau and 
Pomeranchuk [31] . In Subsection 13.31 we will see how to deal with the ultraviolet cut-off by 
invoking renormalization arguments from Quantum Electrodynamics. 

Concerning the description of the electronic structure with A^ electrons, recall that they are 
described by the minimizers Qn of the reduced Bogoliubov-Dirac-Fock energy S^BDF * n the 
charge sectors 

Q A (N) = {Qe Q A , s.t. tr p - (Q) = N}. 

m,0 

Solving this minimization problem is more involved than the construction of the polarized vac- 
uum. The main difficulty arises in the fact that the charge sectors Q\(N) are not stable under 
weak convergence in (3 A . To our knowledge, the construction of the electronic structure with 
N electrons remains open for the Bogoliubov-Dirac-Fock model with exchange term. However, 
Hainzl, Lewin and Sere [22] proved the following Hunziker-van Winter- Zhislin condition (for both 
the reduced and original models). 

Proposition ([22]). Let A > 0, a > and v £ C(R 3 ), and set 

E^bMq) = n inf ^bdf(Q), (23) 

for any real number q. 

(i) Let gGl. // 

Kbdf (?) < Kbdf (?-*) + £°bdf (*) , (24) 

for any fcGl*, then the minimization problem (|23[) owns a minimizer Q q . 

(ii) In case of existence, a minimizer Q q is solution to the self- consistent equation 

lq = Qq + P~ >0 = X(-oo, M9 ) {D q ) + S q , (25) 

where [i q £ [—1,1] is the Fermi level of the electronic structure. The Fock operator D q is defined 
as 

D q = A»,o + a(p Qq - v) * — . 

The self-adjoint operator 5 q is finite rank when fi q 7^ ±1, trace class otherwise, and it ranges in 
the kernel of the operator D q — [i q I . 



When q = N is a positive integer, the polarized vacuum is interpreted as the spectral projector 
l]f C = X(-oo,o](Dn), while the electronic structure with A^ physical electrons is described by the 
spectral projector 7^ = X(o,n N ] (-^Jv)- When the electrostatic potential induced by the charge 
density v is weak enough, that is under condition (|20p . the equality 

tr ^,o (^ C - P -,o)=0, 

holds, which means that the polarized vacuum is neutral. In this case, 

tr P - Q K h ) = tr p - q ( 7Ar - P- ) = tr p - o (Q N ) = N. 

In other words, the electronic structure really exhibits A^ physical electrons which fill in positive 
energy levels of the Fock operator D^. 
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In contrast, nothing rules out the possibility that q is a negative integer. In this situation, the 
minimizer Q q describes a positronic structure with \q\ positrons. The Fermi level \x q is negative. 
The polarized vacuum is identified as before to the spectral projector 7™ c = X(-oo,o](Aj)j while 
the \q\ positrons are represented by the spectral projector 7q h = X(n q ,o]{D q ). 

Checking the validity of condition (j24|) for the Bogoliubov-Dirac-Fock model with exchange 
term is widely open except for some weak-coupling or non-relativistic regimes (see [22] for more 
details). As a consequence, the existence for a given value of N of electronic structures with N 
electrons also remains an open problem. For the reduced model, it is possible to characterize 
precisely the numbers for which inequality (|24p is fulfilled. 

Theorem 1 (\TT\). Let A > 0, a > and v G L 1 ^ 3 ) n C(R 3 ). Set 

Z= [ v. 

(i) There exist two constants q m E [— oo,+oo] and qm G [g m ,+oo] such that the minimization 
problem (I23p owns a minimizer Q q if and only if 

q m <q< qm- 

(ii) Let 

q* — trp— (Q*), 

m,0 

where is the absolute minimizer of the energy <?^b DF on Qa. Then, 

(q*,Z) G [q m ,qM] 2 - 

In particular, when condition (|20p is fulfilled, the minimization problem (|23p owns a minimizer 
Q q for any value of q G [0, Z\. 

When the electrostatic potential induced by the charge density v is weak enough, Theorem Q] 
guarantees the existence of electronic structures with N electrons for any integer N between and 
Z. This is physically relevant in the sense that structures with iV electrons are experimentally 
observed for N between and Z + 1. On the other hand, Theorem [T] does not rule out the 
existence of electronic structures with an arbitrary number of electrons. To our knowledge, there 
are indeed no available upper bounds on \q m \ and qu- 

Computing such ionization bounds is quite involved even for classical Hartree-Fock models 
(see |46| 147] ). Concerning the reduced Bogoliubov-Dirac-Fock model, an additional difficulty lies 
in the sharpness of the ultraviolet cut-off A. Imposing such a sharp cut-off amounts to replacing 
the free Dirac operator D m $ by the operator with Fourier transform 

D^ Q (p) = (a-p + mP)(l + x(|p| 2 /A 2 )) , (26) 

where 

x( r ) = if < r < 1, and x( r ) = +°° if J* > 1. 

The discontinuity in the Fourier transform of is one of the source of troubles which prevents 
from computing upper bounds on \q m \ and qM- 

A natural way to by-pass the difficulty consists in replacing the function \ i R (|26p by a smooth 
function. This has no major consequences on the previous analysis of the reduced Bogoliubov- 
Dirac-Fock model (see |17] for more details). In particular, Theorem [T] remains available with a 
smooth ultraviolet cut-off. Moreover, for the special choice 

X(r) = r 2 , 
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and for a, a\\v\\c and a In A small enough, it is possible to compute the bounds 

- K a ^ A <q m <0, (27) 

and 

Z < q M < 2Z + K a , U:A . (28) 
At the non-relativistic limit a — > 0, A — > +oo and a In A — > 0, the positive number K a ^ Ut \ satisfies 

K a ,v,A -> 0. 

In this limit, one recovers the bound computed by Lieb |35] in the classical case, that is 

q m = and Z < qM < 2Z. 

We refer to [I7| for more precise statements about this topic, as well as for the proofs of (|27p 
and (|28p (which are essentially based on the arguments developed by Lieb in [35J). 

3.3 Charge renormalization for the polarized vacuum 

A crucial ingredient in the proof of Theorem [1] is the following proposition. 

Theorem 2 Q17J). Let A > and a > 0. Consider a function v € L 1 (M 3 ) n C(M 3 ), with 

v = Z G M, 

and denote by q m and qM, the two extremal values in Theorem^ for the interval of existence of a 
minimizer Q q of the energy E^^i'D'F ®^ the charge sector Qa('z)- Given any number q £ [Qrm QM~\> 
the charge density pQ q is an integrable function on R 3 . Its integral is given by 

("-«■) = ifi^- (29) 



where 



A__ 2 



Bl = - V 1 H*" (30) 

vr ./o 1 - z 



In view of Theorem |2j a minimizer Q g is not a trace-class operator, at least when q ^ Z. 
Otherwise, the integral of its density would be equal to 



/ 



PQ q = tr Q q = Q, 



which contradicts the fact that the number in (|30p is positive. 

The fact that the minimizers Q q are not trace class generates a difficulty in their physical 
interpretation. The total electrostatic potential Vph which is induced by the external charge 
density v and the electrons represented by Q q , is defined by the Coulomb formula as 

ir / \ f v iv) - pQ q {y) , rQ -n 

V ph (x) = a / : 2 dy. (31) 

Jr3 \x - y\ 

When |x| is large enough, an approximation for the potential Vph is provided by the expression 

V ph (x) ~ A / {v-p Qq ). 

\ x \ JRJ 
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In view of (1291). it follows that 



V * h{x) ~ l + aBl \x\ * a ~W (32) 

when |x| —> +00. Whereas the minimizer Q q is supposed to represent an electronic structure 
with q electrons, the potential V p h does not match the Coulomb formula for a potential induced 
by a total charge equal to Z — q. 

At this stage, one could argue that it is sufficient to take the limit A — > +00 to solve the 
difficulty. This is not the case. The constant is logarithmically divergent when A — > +00. 
One can check that 

B X = SMA)-rtM(2) +A Oji). (33) 

As a consequence, the potential V p h vanishes in the limit A — > +00. This is another sign of the 
collapse of the model in this limit. 

This logarithmic divergence also appears in Quantum Electrodynamics. This difficulty is 
solved by introducing a charge renormalization. Roughly speaking, charge renormalization con- 
sists in accepting the idea that the bare fine-structure constant a in the model is not the fine- 
structure constant a p h which is experimentally observed. The physical fine-structure constant is 
defined according to the formula 



aph = TT^M' (34) 



a 

so that formula fj32 1) matches with the limit at infinity of a Coulomb potential induced by a total 
charge equal to Z — q. 

Modifying the definition of the fine-structure constant affects in turn the definition of the 
potential Vph in the sense that this potential must be equal to 

V ph (x) = a I ^rdy, (35) 
7r3 \x - y\ 

where p p h refers to the total charge density which is experimentally observed. In view of 031 p . 
one has to admit that the value of p p h is equal to 

«ph/Oph = a P h iy ~ PQg)- (36) 



A natural question is to ask for the physical relevance of the quantities a p ^ and p p ^. Using 
Quantum Electrodynamics, one can compute a power series of p p h with respect to a p h (tending 
to 0), and check that the resulting computations are in agreement with physical measurements. 
In the case of the reduced Bogoliubov-Dirac-Fock model, one can ask for a similar property: 
does it remain possible to compute an expansion of p p ^ when a p h —> 0, and to verify the 
consistence of the expansion with respect to the one provided by Quantum Electrodynamics ? 
The answer is positive provided one introduces a multiplicative renormalization as in Quantum 
Electrodynamics . 

As a matter of fact, our model still contains an ultraviolet cut-off A whose value is unknown. 
In Quantum Electrodynamics, multiplicative renormalization consists in fixing the value of A so 
that the replacement of the bare fine-structure constant a by the physical one a p h amounts to a 
change of physical units (see [10J for more details). In practice, the value of a p h is set equal to 

" P h = Z 3 a, (37) 
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where Z 3 is a fixed positive number. The computations of the power series of p p h with respect to 
a p h are made for Z3 fixed. This amounts to considering A as a function of a p h and Z3 according 
to the identity 

a ph B° A = l-Z 3 . (38) 

The coefficients of the resulting power series are surprisingly independent of the value of Z3. In 
other words, the perturbative computation of the density p p h is independent of the ultraviolet 
cut-off A provided it is fixed according to (|38p . Describing perturbatively the polarized vacuum 
or an electronic structure with N electrons does not really require to set the value of A. 

This property remains true for the polarized vacuum when it is described by the reduced 
Bogoliubov-Dirac-Fock model. 

Theorem 3 ([18]). Let m G N and v G L 2 (M 3 ) n C(R 3 ) such that 

[ ln(l + \k\) 2m+2 \v(k)\ 2 dk < +00. 

We denote by p p h(ctph, the physical density associated to the minimizer of the energy 
^rBDF on 2a according to formulae (|34p , (|36p and ()37p 

[i) Let < e < 1/2. There exist two positive numbers K and a, depending only on m, e and v , 
and a family of functions (f n )o<n<m i n L 2 (R 3 ) nC(M 3 ) ; such that 



/W«ph, 30 - J2 ^«ph L2nc < Ka™+\ (39) 



n=0 



for any < a p h < a and any e < Z3 < 1 — e. 

(ii) The function v$ is identically equal to v , while the functions v n are inductively defined as 

n 

vi=U(vq), and ^ n = ^(^n-i) + ^ Fj (u ni , . . . , v n .) , (40) 

j=3 niH \-rij=n—j 

for n > 2. In this expression, U refers to the Uehling operator defined as the Fourier multiplier 
corresponding to the function 



(41) 



drj. 



U(k)-W C Z ^ ^- 12 " 5|fc|2 | Vg±g (|fcl 2 2)ln(^ ITW + 
TTie nonlinear maps Fj(fj,±, . . . , /i,) are egwa/ to £/ie charge densities of the operators 

QjOi,... = J- / n — V ~ II * n n — "T - 

{Hi) In particular, the functions v n do not depend on Z%, but only on the external charge density 
v. 

Beyond the fact that they do not depend on Z%, the values of the densities v n are consistent 
with the perturbative computations of Quantum Electrodynamics. The function vq is equal to 
the external charge density v. This is exactly the total charge density of the system in the non- 
relativistic case. The function u\ induces a Coulomb potential equal to the Uehling potential 
(see 031 09]) given by 

vw*> = ^ . ^ = f / + V -D»(£ + ?)( I ^-^y) * 
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This potential is the first correction of the polarized vacuum density which is computed by 
Quantum Electrodynamics. 

The proof of Theorem [3] is based on equation (|19p . Assuming that condition (|20p is satisfied, 
the operator <5* identically vanishes, so that we can invoke the Cauchy formula to write 

Here, the operators Q A j are given by 

Q^j = -Tr t> — -— n A(^-PQ») * r - r n A — — — drj. 

2vr y_ oc L> mi0 + zr? V "* | • | D mfi + 177/ 

Their dependence with respect to the ultraviolet cut-off A is explicit through the truncation 
operators Ha defined as 

Ha(/) = /1b(o,A)> 

for any / £ L 2 (M 3 ,C 4 ). 

Translated in terms of the Fourier transforms of the densities pQ t and pq a ., expansion (|42p 
writes as 

PqM = -aB A (k)[pQM - 9(h)) + F A (a(u - PQt ))(k). (43) 
In this expression, the Fourier multiplier B A is equal to 



7T 



i fZ A (\k\) Z 2 -^ \k\ r z ^\ k \) z-4 



B A(k) = - / — ^- — dz + — - * ... dz 



i_ z 2)( 1 + M!(i_ z 2)) 2ttJ yr^A^_M z 



where we have set 



VTTAl- ^/l + (A-r)2 

z A( r ) = : • 



It is useful to write the function B A as 

B A (k) = B° A -U A (k), 

where B° A = B A (0) is defined in §U§ (see [38]). 
The nonlinear map Fa in (|42p is defined as 

= J2 FA > 2n + 1 ^--->v), ( 44 ) 

n>l 

where i*A,n (/•*!> • • • ;/%) stands for the charge density of the operator 



T Ajn (jj,i, . . . , fi n ) = — / — ■ — TT ( IIa * i r -^a ~p: . 

2vr J_ 00 L> mi0 + irj f± V | • | D m)0 + zt/ 



In particular, the functions i*A,n identically vanish when n is even (see |14]). 
In the physical variables a p h and p p h, equation (|42p reduces to 

(l - a ph U A (k)) p^ h (k) + -F A (aphPph) (fc) = ?a(A;), (45) 

with is\(k) = v(k) 1_b(o,2A)(^)- At this stage, it is possible to substitute in (|45p the function p p ^ 
by the formal expansion 

Pph = Yl "A.nQ^h, (46) 
n>0 
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and compute the value of the coefficients VAn- They are inductively given by 

n 

VA,0 = Vh, v\,i = Wa(^a), and v A ,n = ^.(^A.n-l) + F\j(vA ini , . . . ,VA, nj ), 

j=3 niH \-rij=n—j 

(47) 

for any n > 2. In this expression, Z^a refers to the Fourier multiplier associated to the function 
U A - 

Expansion (|39p is then proved into two steps. The formal expansion in (|46p is first rigorously 
derived as a Taylor expansion of order m. This is summarized by the inequality 

m 

Pph - £ ^naShf < Ka™+\ (48) 

n=0 

for a p h small enough. On the other hand, one can check that the coefficients va n converge at 
the limit A — > +oo to the functions v n defined in (|40p . This follows from the convergences of 
the functions Ua and Faj to U, respectively Fj, in the limit A — > +oo. In particular, one can 
replace the coefficients i>A,n in f|48p by the functions v n , so as to obtain expansion (|39|) . 

However, we are not interested in the limit A — > +oo, but in the limit a p h — > 0, with an 
ultraviolet cut-off A fixed so that Z3 = 1 — a p h-BA remains constant. This last assumption is the 
crucial ingredient in the proof. In view of (|33p . this amounts to assuming that 



A^exp^, (49) 

when Qfph — > 0. As a consequence, the Taylor series with respect to a p h of any negative powers 
of A identically vanish. In other words, terms controlled by inverse powers of A play no role in 
the expansion of p p h with respect to a p h. In particular, one can check that 

\Wa,ti ~ Vn\\L 2 nC < Ka ph ln ^ 

for any < n < m. Expansion (|39p follows combining with (j48]). 

To conclude the derivation of Theorem [3l notice that the dependence on Z% is entirely con- 
tained in 049 p . so that it vanishes when a Taylor expansion of an inverse power of A is performed. 
This explains why the coefficients v n do not depend on Z^. 



4 The Pauli-Villars regulated model 
4.1 Formal derivation 

A major restriction in the (reduced) Bogoliubov-Dirac-Fock model lies in its purely electrostatic 
nature which prevents, for instance, from describing the role played by photons in the vacuum 
polarization. Our goal is now to explain, at least formally, how to derive a more general Hartree- 
Fock model taking into account some features of the photons as well as of external magnetic 
fields. We refer to |16j for more details (see also [33J). 

The main difficulty arises in the choice of the regularization which we have to introduce in 
order to define the model properly. Recall that the formal Hartree-Fock energy for describing 
the polarized vacuum in the Coulomb gauge may be written according to (|5]) as 

*m=i.(wj(7-v>))4/ / 

v / z j K 3 j R 3 \x — y\ 

e 2 [ [ |(7-l/2)(a;,y)| 2 , , If, , „, Nl2 , 

/ / — — , 1 — 1 dxdy + — / \c\\AA{x)\ 2 dx. 

2 Jr3 J R 3 \x - y\ 8vr J R3 
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The derivation of the (reduced) Bogoliubov-Dirac-Fock model from this energy consists in omit- 
ting some terms on one hand, introducing a regularization on the other hand. The regularization 
is based on the introduction of an ultraviolet cut-off A in the Fourier space. This choice breaks 
the magnetic gauge invariance corresponding to replace A by A + Vy?, which is used in Quantum 
Electrodynamics to eliminate some divergences in the perturbative computations. As a result, a 
relevant model including photons and external magnetic fields cannot rely on a sharp ultraviolet 
cut-off. 

In order to derive an alternative model, it is convenient to express the direct electrostatic 
term in ([5]) in terms of the Coulomb potential 



\x - y\ 



according to the formula 



e 



p 7 _ 1/2 (x)p 7 _ 1/2 (y) f 1 



dxdy = e p 7 _i/ 2 (^) K-1/2O) dx - — / \VVy_ 1/2 (x)\ dx. 
\x-y\ i R a 1 11 8vr7 R 3 ' 

Since the potential solves the maximization problem 

el p y -i/2(x)V r - 1 / 2 (x)dx- — I |Vy 7 _ 1/2 (x)| 2 (ix = sup je / p 7 -i/ 2 V - — [ |W 

we can write the formal energy in ([5]) as 

£^( 7 ,A) = sup ££r(7,A). 
y 

where we have set A = (V, .A). In this formula, the formal Hartree-Fock Lagrangian C^fp is 
given by 

C^,A)=tr(D m>e{A+Ae ^-l/2))-^ [ I 1(7 "i 1 / 2 ^'^' 2 dxdy 

z Jm 3 Jr 3 \ x — y\ 

+ — f (|curlA(x)| 2 - |VV(x)| 2 ) dx. 
8tt Jr3 

We can also introduce a reduced Lagrangian omitting the exchange electrostatic term according 
to the formula 

^(7,A)=tr(£» m)e(A+Aext) ( 7 -l/2)) + ^- f {\cuilA(x)\ 2 -\VV(x)\ 2 )dx. 

87r Jr 3 

In the sequel, we restrict our attention to the reduced formalism. To our knowledge, the analysis 
of the original model remains largely open (see |16j for more details). 

The polarized vacuum is constructed as a minimizer of the formal reduced energy 

Sfe(7,^)=sup £^(7, A), 
v 

with respect to the one-body density matrix 7 and the classical photon field A. The operator 
7 is an orthogonal projector as before, while the field A is divergence free due to the Coulomb 
gauge. In case of existence, a minimizer (7*, A*) is solution to the self-consistent equations 

7* = X(-oo,0](An,e(>l,+Aext))> 

-AA* = 47rej > _ 1 / 2 , 
-AT4 = 47rep 7>fc _i/ 2 , 
div A* = div A ext = 0. 
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In this expression, the charge density p 7<1 _i/2 an d the charge current j 7it _i/2 are given by 

P-y,-i/2{x) = tr C 4 ((7* - \/2){x,x)) and j 7 ,_i/ 2 (x) = tr C 4 (0(7* - l/2)(x,x)), 

where (7* — 1/2) (x,y) refers as above to the (formal) kernel of the operator 7* — 1/2. The 
first equation in (|50p is consistent with the picture of the Dirac sea since the minimizer 7* 
is the negative spectral projector of the Fock operator -D m . e (A»+A cxt )- The equations for the 
electromagnetic four-potential A* = (V*,A*) are well-known in the Physics literature (see e.g. 

mi)- 

Regarding the construction of electronic structures with iV electrons, a charge constraint of 
the form 

tr (7 - 1/2) = N, 

is added as for the reduced Bogoliubov-Dirac-Fock model. The equations for the minimizers 
(77V, An) write as 

77V = X{-oo,ti N ]{D m ,e{A N +A cx t))> 

-AA N = 47rej 7Ar _ 1/2 , 
-AV N = 47re/) 7jv _i /2 , 
div A N = div^ ext = 0, 

where /i/v is the Fermi level of the electronic structure. According to the picture of the Dirac 
sea, the physical electrons are identified to the projector 7^ h = X(o,fj, N ](D mt e(A N +A cxt ))i while 
the virtual electrons of the polarized vacuum are represented by 7™° = X(-oo,o] (^m,e(Ajv+A cxt ))- 
The model describes a positronic structure with \N\ positrons when N is a negative integer. In 
the sequel, we focus on the construction of the polarized vacuum. Constructing electronic or 
positronic structures also remains an open problem. 

Instead of maximizing the Lagrangian £^^,(7, A) with respect to V and then minimizing 
the resulting quantity with respect to 7 and A, one can follow the alternative strategy which 
consists in minimizing first with respect to 7, and then looking for a saddle point with respect 
to V and A. We do not claim that we are solving exactly the same problem in this way, but 
this alternative strategy provides a polarized vacuum which is consistent with the self-consistent 
equations (|50p (see Theorem 2]). Moreover, the problem is simplified. The unique minimizer of 
the minimization problem 

inf£^( 7 ,A), (51) 

is indeed explicitly given by 

7A = X(-oo,0](An,e(A+A cxt )), 
with a value for the minimum equal to 

inf £^( 7) A) = -\tr\D mAA+Acxt) \ (I curlA| 2 - |W| 2 ). 

In order to construct the polarized vacuum, it only remains to solve a min-max problem which 
only depends on V and A, namely 

nun max { - i tr \D mAA+Acxt) \ + -L J (| curM| 2 - |W| 2 )}, 

where |T| stands for the absolute value of the operator T. 

At this stage, it is necessary to acknowledge that most of the previous discussion is only 
formal. We now have to introduce regularization techniques from Quantum Electrodynamics 
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in order to define a rigorous model. The first element in this direction is reminiscent from the 
(reduced) Bogoliubov-Dirac-Fock model. When A = A cx t = 0, the negative spectral projector 
P~ of the free Dirac operator formally solves the problem (|5ip . However, the value of the 
minimum is equal to 

inf £° HF (7,0) = --tr |A»,o|, 
7 2 1 1 

which is an infinite quantity. In order to deal with finite quantities, and since the situation under 
our consideration corresponds to the free vacuum, we define a relative energy according to the 
formula 

^r(A) = ~ tr (|An,o| - \D mMA+Acxt) \) +^j R3 {\ curlvll 2 - |W| 2 ), 

so that the free vacuum now has a vanishing energy. Since this amounts to adding an infinite 
constant, this does not change the variational problem under our analysis, but we can now hope 
that the quantity f^f™* (A) is finite provided A and A ext belong to some suitable function space. 

Actually, the energy remains divergent for large momenta. As a matter of fact, the operator 
l-^m.ol — l^m,e(A+A cxt )l i s n °t trace class when A + A ext ^ 0, so that its trace is not well- 
defined (see [37] for more details). In order to remove these divergences, an ultraviolet cut-off 
is necessary. Various techniques from Quantum Electrodynamics are available. Our choice is to 
rely on a regularization proposed by Pauli and Villars in [39J, which consists in introducing the 
functional 

£^{A) = ±tr \J2°S {\ D mJ ~ \D mj , e(A+AeKt) \fj +^J Ki (|curlA| 2 - |W| 2 ). (52) 

In this expression, the index j = corresponds to the physical electron-positron field. In partic- 
ular, mo = m is the (bare) mass of the electron. The indices j = 1 and j = 2 describe fictitious 
heavy particle fields. Their role is to remove the worst ultraviolet divergences (the linear ones 
to be more precise). In order to reach this goal, it is well-known (see [39J) that the coefficients 
Cj and the masses rrij must fulfil the two conditions 

J J 
j=o 3=0 

At least two additional distinct masses m\ and 777,2 are therefore necessary. In the sequel, we fix 
J = 2 and cq = 1. For this choice, the condition (|53p is equivalent to 



777n — ITlo , 7777 — 777 rj 

ci = — 5 5- and C2 = — 5 

7772 — m l m 2 ~ m l 

In the limit m\ — > +00 and 7772 —> +00, the regularization does not prevent a logarithmic 
divergence which is identical to the divergence of the constant B\ in (|33p (see Proposition [1]) . 
The divergence is best understood in terms of the averaged ultraviolet cut-off A defined as 

2 

log(A 2 ) = -^ Ci log(777 2 ). (54) 
3=0 

The value of A does not uniquely determine the masses mi and 7772. In practice, they are chosen 
as functions of A such that the coefficients c\ and C2 remain bounded when A — > +00. 
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In the Coulomb gauge, the Pauli-Villars regulated energy <?py xt is rigorously well-defined under 
the natural conditions that the fields B = curlA, B cxt = curl A ext , E = — W and E cxt = — W cx t 
are square integrable (see Proposition [I]) . The polarized vacuum is described using a solution to 
the min-max problem 

minmax£'py xt (A). (55) 
More precisely, it is identified to the projector 

7PV = X(-oo,0] (- D m ,e(A,+Aext)) 

where A* = (V*,A*) is a saddle point for the problem (|55p . One can guarantee the existence of 
such a saddle point at least when the external electromagnetic field _F cxt = (-E7 ex t ; B cxt ) is weak 
enough. 

4.2 Construction of the polarized vacuum 

We now define properly the Pauli-Villars regulated energy in (|52p before solving the min-max 
problem (|55p . The natural framework for defining the energy Spy* 1 is provided by the Coulomb- 
gauge homogeneous Sobolev space 

#] iv (R 3 ) = |A = (V,A) G L 6 (R 3 ,R 4 ), s.t. div A = and F = (-W,curlA) G L 2 (R 3 ,R 6 )}, 

which is an Hilbert space for the norm 

II^IIjL (rs) = II w IIl 2 (r 3 ) + II curl ^lll 2 (R3) = \\ F \\h 



When A G -ffj iv (R 3 ), the integral in (|52p is well-defined, but we have to provide a rigorous 
meaning to the first term in 052 p . This amounts to defining properly the functional 



, 2 

T PY (A) = - tr ^2 C J ( | D mj ,o | - | D mj tA \\ (56) 



2 

tr" 

2 

for an arbitrary four-potential A G i/j iv (R 3 ). In this direction, we can establish the next 
proposition. 

Proposition 1 (|16j). Assume that the coefficients Cj and the masses m,j satisfy 

2 2 

Co = 1, m% > mi > rriQ > and Cj = c j m2 j = 0- (57) 

j=0 j=0 

(i) Let 

H 2 

T A= g £^(1^,01 -\D mj ,A\). (58) 

j=0 

1 ("in>3\ ill nn>3\ j.u„ +v _ . T . „•„ ™ r2 



Given any A G L (R )f"l.fZjj v (R ), the operator tr 'jm T a is trace class on L (R ,C). In particular, 
the quantity .Fpv(A) is well-defined by the expression 

J-pv(A) =tr(tr C 4T A ). (59) 



n) T/ie functional J-py can be uniquely extended to a continuous mapping on H^ iv ( 

T PV (A)=T 2 (F)+K(A), (60) 



(Hi) Let A G #] iv (R 3 ). We have 
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where F = (E, B), with E = — W and B = curl A. The functional TZ is continuous on H^iv 
and satisfies 

\ v j=0 v 3=0 3 / 

/or a universal positive number K. 

(iv) The functional Ti is the non-negative and bounded quadratic form on L 2 (M 3 ,M 4 ) given by 



where 



MF) = ^J Ra M(k)(\B(k)\ 2 - \E(k)\ 2 ) dk, 
2 2 rl 

M(k) = ^2 cj / u(l - u) log (m 2 + u(l - u)\k\ 2 ) du. 



(62) 



(63) 



The function M is positive and satisfies the uniform estimate 



where A is given by ([54 



< M{k) < M(0) = 21 ° o g(A \ (64) 

67T 



The proof of Proposition [I] relies on a perturbative expansion of the operator Ta- Invoking 
the formula 

T = — / 2 - — — + — — <kj, 
2tt 1 + iuj 1 — iuj J 

the operator Ta may be written as 

1 / \ / 2W \ 

Ta = — / > c 7 ' ( 1 duj. 

47r JM. V ^m, , A + iU Dm, ,A ~ i(*> D mj >0 + IUJ D mj fl ~ IUJ J 

Expanding with respect to the powers of A leads to the expression 
5 

T A =Y, T n( A )+n{A) 

n=l 

=J-E / {Rn(u,A) + R n (-oj,A))doj + ±- I (R' 6 (uj,A)+R' 6 (-uj,A))duj, 

with 

2 

IUJ 

3=0 



it! n (a;,A) = ]>> ~_. ((tt-A-y) - 1 , 

for 1 < n < 5, and 



i^(o;,A) = Vcj - l - . ({cc-A-V) - 1 , )*. 



3=0 

'2 



Due to the conditions (|57p . the operators tr C 4 T n (A) and tr C 4Tg(A) are trace class on L v 
when A belongs to L 1 (R 3 ) n ifj iv (M 3 ). The quantity .Fpv(A) in (|59|) is therefore well-defined. 
Moreover, it depends Holder continuously on A £ i/j iv (M 3 ), so that it can be extended to the 
space if i iv (M 3 ). 
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Let us emphasize the introduction of the C 4 -trace here. The operators T n (A) and 7g(A) are 
probably not trace class without taking first the C 4 -trace (except when ^4 = 0). Defining J-pv 
as in (|59p extends the formal definition (|56p to the case where Ta is not a trace-class operator. 
The two definitions remain identical when Ta is trace class. 

Concerning the second order operator tr C 4 T%(A), an explicit computation leads to the formula 

tr(tr c4 T 2 (A)) = F 2 {F), 

where J-2(F) is defined in (|62p . The Fourier multiplier M in (j62]) describes the linear response 
of the virtual electrons in the polarized vacuum. In view of the convergence 

l im _ M(k)) = U(k), (65) 



where U is the Uehling multiplier given by (|4ip . the function M appears as the Pauli-Villars 
equivalent of the function in the reduced Bogoliubov-Dirac-Fock model. This similarity in 
the two models results from the gauge and relativistic invariances in Quantum Electrodynamics. 

It follows from 065 P that the self-consistent equations of the charge and current densities 
corresponding to a (possible) solution A* to the min-max problem (|55p are very similar to the 
equation f|43[) for the charge density pQ t of a minimizer Q* of the reduced Bogoliubov-Dirac- 
Fock model. Even if this was not done in |16j . the renormalization technique applied to define 
a physical charge density p p h and to compute its perturbative expansion with respect to a p h 
in the case of the reduced Bogoliubov-Dirac-Fock model is likely to work the same with the 
Pauli-Villars regulated model corresponding to the energy £p^ xt . 

This however requires to construct first a solution A* to the min-max problem (|55p . The 
construction is possible when the external electromagnetic four-potential A ext is small enough. 
In this case, one can deduce from the expression of the second-order functional T% in (|62p that the 



energy £py xt owns a local saddle point geometry close to the four-potential A = 0. The existence 
of a (local) solution to the min-max problem (|55p follows using tools from convex analysis. 

Theorem 4 ([E]). Assume that the coefficients Cj and the masses rrij satisfy the conditions 

(EZD- 

(i) There exists a positive number r such that, given any four-potential A ex t G -ffj iv (K 3 ) such 

that 

rJmQ 

e||-A«t||ji v (R3) < — — , (66) 
there exists a unique solution A* G i7j iv (M 3 ) to the min-max problem 



(67) 



£^ xt (A*)= max inf £p^{A) 

II V ^ / IIl2<^ J1 \\curlA\\ L2 <^L-0 

= min sup £j^y xt (A). 

||curlA|| i2 <^ || W||i2< M 

(ii) When A ext = 7 the solution A* is equal to 0. 

(Hi) The four-potential A* is a solution to the nonlinear equations 

[~^ v r A r p *> (68) 

where G C(R 3 ) and G C(1R 3 ) are defined as 

p*(x) = [tr C 4 Q*] (x, x) and Ja{x) = [^c A a Q*] (%,%)■ (69) 
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In this expression, the function Q*(x,y) refers to the kernel of the locally trace-class operator 

2 

Q* = ^ c iX(-c»,0](- D m J -,e(A.+i4« rt ))- (70) 
3=0 

According to the previous derivation of the Pauli-Villars regulated energy, the polarized vac- 
uum is identified to the projector 

7PV = X(-oo,0] (^m ,e(A»,+A cx t))- 

Its construction is only local and only available for small enough external electromagnetic fields. 

To our knowledge, the existence of a global solution to the min-max problem (|55p remains an 
open problem. A first attempt to answer this question could concern the property that A = is 
the unique global saddle point of the energy Spy. 

The construction of a (local) minimizer in large external electromagnetic fields is another 
appealing problem, in particular since it certainly requires to understand the phenomenon of 
production of electron-positron pairs (see |42j for a first analysis of this phenomenon). 
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